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Abstract 

Recent developments of dynamic x-ray characterization experiments of dense matter are re- 
viewed, with particular emphasis on conditions relevant to interiors of terrestrial and gas giant 
planets. These studies include characterization of compressed states of matter in light elements by 
x-ray scattering and imaging of shocked iron by radiography. Several applications of this work are 
examined. These include the structure of massive "Super Earth" terrestrial planets around other 
stars, the 40 known extrasolar gas giants with measured masses and radii, and Jupiter itself, which 
serves as the benchmark for giant planets. 
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I. INTRODUCTION 



We are now in an era of dramatic improvement in our knowledge of the physics of mate- 
rials at high density. For light elements, this theoretical and experimental work has many 
applications, including internal confinement fusion as well as the interiors of gas giant plan- 
ets. For heavy elements, experiments on silicates and iron at high pressure are helping to 
better understand the Earth, as well as terrestrial planets as a class of objects. In particu- 
lar, the discovery of rocky and gaseous planets in other planetary systems has opened our 
imaginations to planets not found in our own solar system 

While the fields of experiments of matter at high densities, first principles calculations 
of equations of state (EOS), planetary science, and astronomy do progress independently of 
each other, it is important for there to be communication between fields. For instance, in 
the realm of planets, physicists can learn of key problems that exist in the area of planetary 
structure, and how advances in our understanding of input physics could shed new light in 
this area. Astronomers and planetary scientists can learn where breakthroughs in physics 
of materials under extreme conditions are occurring, and be ready to apply these findings 
within their fields. 

This brief review focuses on work presented at the joint American Physical Society (APS), 
High Energy Density Laboratory Astrophysics (HEDLA), and High Energy Density Physics 
(HEDP) meeting in April, 2008. We first discuss some experimental and theoretical work on 
light elements, including some applications to gas giant planets, which are predominantly 
composed of hydrogen and helium. We discuss new models of the interior structure of Jupiter 
and review the observed mass-radius relationship of extrasolar giant planets (EGPs). We 
then focus on terrestrial planets and investigate new experimental work on iron. We then 
discuss the predicted structure of "Super Earth" planets, predominantly solid planets made 
up of iron, rock, and water, from 1-10 Earth masses (M^). 

II. LIGHT ELEMENTS AND GIANT PLANETS 
A. Pulsed x-ray probing of light elements 

For accurate measurements of densities and temperatures in dense and compressed mat- 
ter, novel x-ray Thomson scattering techniques^] have been developed. These experiments 
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employ powerful laser-produced x-ray sources that penetrate through dense and compressed 
materials with densities of solid and above. Both thermal Ly-a and He-a radiation from 
nanosecond laser plasmas 3] or ultra short pulse laser-produced K-a inner-shell emission j4| 
have been shown to fulfill the stringent requirements on photon numbers and bandwidth 
for spectrally resolved x-ray scattering measurements in single shot experiments. Experi- 
ments have been performed in the non-collective (backscatter) regime, where the scattering 
spectrum yields the Compton feature j^. On the other hand, collective scattering on plas- 
mons, i.e., electron density (Langmuir) oscillations, have been observed in forward scatter 
geometry p]. 

These techniques have recently been applied to shock-compressed beryllium (H. J. Lee 
et al., in prep). In experiments at the Omega laser facility 0] twelve laser beams (500 J 
each) directly illuminate the foil with laser intensities of J = 10^^ — 10^^ W cm~^ producing 
pressures in the range of 20-60 Mbar and compressing the foil by a factor of 3. The Compton 
scattering spectrum of the 6.18 keV Mn He-a and 6.15 keV intercombination x-ray probe 
lines measured at G = 90" scattering angle shows a parabolic spectrum downshifted in energy 
from the incident radiation by the Compton effect; the shift is determined by the Compton 
energy Ec = h^k^/47rme = 74 eV, with k = 47r(Eo/hc) sin(6/2) = 4.4A"\ and Eq the 
energy of the incident x-rays. The Compton scattering spectrum directly reflects the electron 
distribution function; for a Fermi-degenerate system the width of the Compton spectrum 
provides the Fermi energy, Ep ~ n^/^. Unlike for plasmas with a Maxwell-Boltzmann 
distribution, the width is sensitive to the electron density. In addition, the intensity ratio of 
the elastic to inelastic scattering feature from Fermi-degenerate plasmas is sensitive to the 
ion temperature because elastic scattering is dependent on the ion-ion structure factor. 

Figure ([1]) shows the scattering data along with calculated scattering spectra for which 
the electron density (left) and the temperature (right) has been varied. For the analysis 
we assume Te=Tj and Z = 2 consistent with calculations and with the measurements from 
isochorically heated Be. Density and temperature obtained in this way are ne = 7.5 x 10^^ 
cm~^ and T = 13 eV for representing a Fermi temperature of Ep = 30 eV and scattering 
parameter a = l/{kXs) = 0.48. The error bar for the measurement is of order < 10%, 
dominated by noise. This error estimate is not affected by uncertainties of Z because the 
shape of the Compton scattering proflle provides an additional constraints. The parameters 
inferred from the theoretical flt match radiation hydrodynamic simulations of this experi- 
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ment to 10% and also agree with the results of forward scattering measurements that have 
independently measured Ue and Tg from the plasmon spectrum. 

The experiments have directly measured the conditions and dynamic structure factors 
of shocked matter, thus going beyond characterization of shock wave experiments with par- 
ticle and shock velocities. These novel experiments have only now become possible with 
the advent of penetrating powerful x-ray probes produced on high-energy density physics 
facilities. This feature further allows testing of radiation-hydrodynamic calculations with 
different EOS models for shock- compressed matter. Future experiments will apply Comp- 
ton scattering to measure the compressibility and adiabat of compressed matter, including 
hydrogen jsl. 



B. Experiment and Theory of H/He 

Experiments on deuterium and helium provide vital EOS data to model the interior of 
giant planets (Fig. [2]). Recent work include explosively driven shocks in hydrogen and 



deuterium 



M] that provided experimental evidence of pressure-driven dissociation 



transition 

in dense hydrogen that is associated with a rapid increase in conductivity 0, Q. The 
results in Fortov et al. [9] were interpreted as signs of a plasma phase transition of first- 
order. However, within experimental error bars, an interpretation in terms of a gradual 
dissociation transition as predicted from first-principles simulations [11] is also possible. 
While standard shock wave experiments can reach very high temperatures and megabar 



pressures, they increase the sample density only 4 to 5 fold [12|. There is consequently 
very little data taken directly under the conditions found inside Jupiter and Saturn. In 
addition to existing isentropic and off-hugoniot shock compression techniques, this issue is 
addressed with a new experimental method that combines static and dynamic compression 



techniques. By precompressing the sample statically in a modified diamond anvil cell 13| . 



a higher final shock density is reached, and one can therefore probe deeper in the planetary 
interiors. In Figure ([3]) we show a comparison of the Eggert et al. |3] precompressed data 
(the hashed region in Fig. [2]) to first principles simulations of Militzer [ij] . The agreement 
is particularly good at high compression but deviations are observed for the measurements 
without precompressions. Further experimental and theoretical work will be needed to 
resolve this discrepancy. 
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C. Jupiter and Extrasolar Giant Planets 



An important requirement of theories of planet formation is to account for the present 
day structure of Jupiter and Saturn, which are relatively well observed. The history of their 
formation is imprinted primarily in the amount and distribution of heavy elements in their 
interior. Heavy elements are supplied as solid bodies while the massive H/He envelopes of 
giant planets accumulate through gas accretion. The relative importance of the two accretion 
processes during the formation of the planet is closely tied to the formation process and the 
surrounding environment. As noted by Saumon and Guillot [l5|, uncertainties in H/He EOS 
dominate all uncertainties when trying to understand the interior structure of Jupiter. It is 
then essential that accurate EOS measurements can be made for light elements under giant 
planet conditions. 

Interestingly, the first two modeling efforts based on first principles EOS computed inde- 
pendently but with essentially the same method give very different results for the amount 



and distribution of heavy elements in 
discussed in Militzer and Hubbard 



upiter Uj, ll6|] . The reasons for this discrepancy are 
17( 1 and stem primarily from different assumptions for 
the interior structure of Jupiter. Nettelmann et al. 16j assumed a different concentration of 
heavy elements and helium for the molecular and the metallic regime that could for example 



be introduced by a first order phase transition in hydrogen. Militzer et al. 



ll| found no 



evidence of such sharp transition in their first-principle simulation and concluded the mantle 
must be isentropic, fully convective, and of constant composition. Jupiter's interior struc- 
ture, as derived by Militzer et al. [llj, is shown in Figure (jlj). We may be at the threshold 
where the EOS of H/He mixtures is understood well-enough to force qualitative changes in 
our picture of the basic structure of Jupiter. The implications for the interior of Saturn and, 
by extension, for the planet formation process remain to be explored. 

New observational data on Jupiter is hard to come by, given that space probes are nec- 
essary to measure deep atmospheric abundances and map the gravitational and magnetic 
fields. In 2011 NASA will launch the JUNO orbiter, which will reach Jupiter in 2016 18|. 
This orbiter has several important goals relating to the structure of the planet. The deep 
abundance of oxygen will be measured, which is potentially Jupiter's third most abundant 
element, after H and He. The detailed mapping of Jupiter's gravity field will give us unpar- 
alleled access into the internal structure of the planet 191 ]. 
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The available data set on transiting EGPs continues to expand. These planets periodically 
pass in front of their parent stars, allowing for a determination of planetary radii. Planetary 
masses are determined from the Doppler shift of the parent star's spectral lines. Forty such 
planets, most in very close-in orbits, are now known around other stars. The doubling 
time for the number of detected planets is now less than one year [20]. Dedicated space 



missions such as CoRoT 



21 



22j and Kepler [23| will detect dozens to hundreds of additional 



Neptune-like to Jupiter-like planets, in additional to smaller terrestrial planets, discussed 
below. In Figure ([5]) we show the measured masses and radii of known planets, compared 
to the solar system's giant planets, and predictions from theoretical models at two different 
irradiation levels [24]. 

It is the extreme diversity shown in Figure that is most surprising. Irradiated giant 
planets were expected to be inflated relative to Jupiter [25], but the range of radii for planets 
of similar masses does not yet have a satisfactory solution. To explain the relatively smaller 
radii of some of the Jovian planets, heavy element abundances of 100-200 (0.31-0.62 
Mj) are needed. However, we are ignorant of whether these heavy elements are predomi- 
nantly mixed into the H/He envelope or within a distinct core. This issue is actually even 
more complicated, due to several additional factors 26| : (i) the differences between the var- 
ious equations of state used to characterize the heavy material (water, rock, iron), (ii) the 
chemical composition of the heavy elements (predominantly water or rock?), and (iii) their 
thermal contribution to the planet evolution (which is often ignored altogether). Deriving 
heavy element abundances based on a given planet's mass and radius will be uncertain, given 
these issues. As we discuss in §111, similar issues are also important for massive terrestrial 
planets. 

To explain the large radii of many of the planets, either an additional internal energy 
source must be invoked, 



27 



or that the cooling and contraction of these planets has been 
St ailed. |29|,l30| The only clear trend to date is that planets around metal-rich parent stars 
tend to possess larger amounts of heavy elements. 31] 
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III. HEAVY ELEMENTS AND TERRESTRIAL PLANETS 



A. Experiments on Iron 

With current technology, diamond cell experiments (static), do not allow one to obtain 
meaningful data at temperatures of several thousands of K once the pressure exceeds 200 



GPa. Although measurements can be achieved with good precision below 200 Gpa [32] , the 
melting curve of iron or iron alloys at the inner core boundary (330 GPa, about 5000 K) is 
thus beyond the capabilities of these experiments. On the other hand, dynamic experiments 
can easily reach inner core pressures but the corresponding temperatures, which are large 
and are fixed by the Hugoniot curve, do not allow one to explore the relevant P-T space. 
As a result, the iron phase diagram at conditions corresponding to the Earth's inner core 
has never been directly measured and large uncertainties remain regarding its equation of 
state (EOS). These unknowns severely limit current Earth modelling as the iron EOS is of 
utmost importance to constrain the chemical composition and energy balance of the Earth's 
core. The discovery of low-mass planets outside the solar system renders the exploration of 
iron at > 500 GPa (5 Mbar) pressures and ~1 eV temperatures even more pertinent. 

The French National Research Agency (ANR) recently funded a several-year program 
focused on the development of new diagnostics to study the physical properties of iron, the 
development of methods to explore broader regions of the EOS diagram, and the combined 
use of experimental and theoretical methods to characterize the high pressure phases of this 
element. 

The first part of this project is to develop adequate x-ray sources both to radiograph 
and/or perform diffraction measurements on shock compressed iron. Results at lower ener- 
gies have already been obtained for aluminum H . To this aim, experiments were performed 
on the 100 TW laser system at LULI, France, which delivers 20 J in 0.3-10 ps at a wave- 
length of 1057 nm and 6J when frequency-doubled. The latter was used to look at the 
effects of preformed plasma due to the laser 'pedestal' which is 500 ps wide with a contrast 
of 10"^. Different target materials and geometries, as well as the effects of laser parameters 
and filtering/shielding of the detector were studied. At very high laser intensities, x-rays are 
generated by energetic electrons produced by the laser-plasma interaction which penetrates 
the target and produces K-a radiation. The x-ray emission stops a few ps after the end of 
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the laser pulse, as the electrons lose their energy due to classical charged-particle stopping 
processes. For laser pulses of ~10 ps, a temporal resolution of less than 20 ps can then be 
achieved. This duration is short enough to resolve Shockwaves for density measurements in 
EOS studies, as the shock velocities are of the order of a few tens fim ns~^. 

The backlighter target was made of W, (producing Ka energy of 60 keV). To produce a 
small source size, required for high 2D spatial resolution, we used thin (18 fim diameter) W 



wires 3J]. The spatial resolution of the x-ray source was measured using a crossed pair of 100 
/im diameter gold wires. The latter was used for these high energy x-rays, as the absorption 
of a standard gold grid was too low (~5%). The magnification of this point-projection 
system was 30 x, with a source-detector (imaging plate filtered with Tm) distance of 30 cm; 
an additional 2.5 cm plastic layer was used to stop energetic charged particles emerging from 
the target. To characterize the spectral distribution of the x-rays, a transmission crystal 



spectrometer (DCS) [35(], was implemented in order to measure the contribution of high- 
energy x-ray background to the radiograph image. The experiment reliably delivered high 
quahty radiographs of static targets with best results obtained at intensities of 10^*^ W cm~^. 
Figure (EI) shows a radiograph at 60 keV of a test target (100 /im thick gold wire at 30 x 
magnification) obtained from a 18 fim diameter W-wire target. Analyzing the absorption 
profile of the wire on the detectors shows a spatial resolution of better than 20 fim. The 
contrast on iron steps (right part of Figure [6]) shows the resolved density gradients that 
allows for the deduction of the density of shocked iron with error bars lower than 10 %. 

This technique has then been applied to radiograph a shock-compressed target made of 
a ablator pusher and a 500 /im diameter iron disk, 250 /im thick. For this experiment, 
besides the short pulse beam, a high energy long pulse beam was needed to drive a uniform 
planar shock. Therefore we used the new LUL12000 facility which has this capability and 
obtained the first radiograph of a laser shock compressed iron target (Figure [7]). Due to 
a lower contrast on the LULI2000 facility than the 100 TW, the signal/noise ratio is not 
optimum but the shock front is clearly observable. Detailed analysis is still underway. 

B. Application: Super Earths — Massive Terrestrial Planets 

Super-Earths are the newest class of discovered extra-solar planets. These 1-10 earth- 
mass (M®) planets are likely to consist of solids and liquids rather than of gases. With their 
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relatively large masses, they experience very large internal pressures. Pressure constrains 
the power law relationship between mass (M) and radius {R) of solid planets. The value for 
the exponent in R = R^ef{M/M^ef)^ is 0.262 < /? < 0.274 as constrained by the different 
internal structure models for super-Earths, while it is /9 = 0.3 for planets between 5-50% 
the mass of Earth [36|, l37|. -Rref is the radius of a planet with reference mass Mref, usually 



Earth's, that may be rocky or have large amounts of H2O with a correspondingly larger 
-Rref- The central pressure of rocky super-Earths (up to ~60 Mbar) scales proportionately 
with mass, reaching values that challenge the understanding of rock behavior under such 
extreme conditions. Despite the different treatments in the models 2S, |38|, l39|] and intrin- 
sic uncertainties in the equation of state (EOS), composition and temperature structure, 
the mass-radius relationship is robust, and thus, useful for inferring the expected signal in 
searches for transiting super Earths. 

However, information on the structure, such as the size and state of the core, crucially 
depends on the exact behavior of super-Earth materials (silicates, iron, iron alloy and ices) 
at high pressures and temperatures. In order to accurately describe the physical properties 
of super-Earths, such as their ability to have a magnetic field by having a molten core, or to 
extract information to constrain formation models such as from the existence of a metallic 
core, we need a very detailed description of super-Earths' interior, that can not be done 
without improvements in the EOS of silicates, iron alloys and ices. 

A few questions that, if addressed, will considerably improve the internal structure models 
and thus, our interpretation of the data are: (1) What is the stability field of post-perovksite 
and are there other higher pressure silicate phases? Our lack of knowledge of other existing 
phases means that the radius in models is an upper value. (2) At the pressure range of 
super-Earths (up to 60 Mbars), which existing EOS is more accurate? A few high pressure 
experiments can illuminate the extrapolation qualities and deficiencies of the different EOS 
used by the models (Vinet, Birch-Murnaghan, ANEOS, Thomas- Fermi-Dirac, etc). However 
significant process has been made with first-principles computer simulation [40| . How much 
iron can post-perovskite accommodate? If this high-pressure silicate phase accommodates 
a large amount of iron (as suggested by Mao et al Mi), it could affect the size of the core 
and to a smaller extent the total radius of the planet. (4) What are the thermodynamic 
properties of all mantle materials, especially the Gruneisen parameter? Post-perovskite has 
a more sensitive Gruneisen parameter to volume than perovskite, such that the temperature 
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profile of a mantle made mostly of post-perovskite would be cooler than that of perovskite. 

In addition, to infer planetary composition from the M and R data that will be available 
in the next few years, we need accurate EOSs. Even without errors in the data and structure 
models, a large number of compositions can fit the same average density (see Figure [H&). 
The uncertainty in radius from EOS is ~2-3%, which will be comparable to the precision 
that powerful space telescopes will yield in follow-up observations (i.e. the James Webb 
Space Telescope, JWST). By reducing the uncertainty in the EOS we make the structure 
models more accurate and useful. Thus, there is a need for accurate equations of state of 
solid planetary materials to pressures up to ~ 60 Mbar, the central pressure of the densest 
and largest (IO-M0) super-Earth. 

IV. CONCLUSIONS 

New experiments are now probing states of dense matter that were previously beyond our 
grasp. In particular x-ray techniques are allowing us a view into materials that previously 
had been hidden. At the same time first principles techniques are allowing accurate 
determination of EOSs for planetary interest. Off-Hugoniot experiments of H and He will 
test these EOS and lead to more accurate models for Jupiter and Saturn. In addition, new 
planets are being discovered at an accelerating rate, which will continue to expand the 
limits of P-T space that are of "planetary" interest. Our current era is one of dramatic 
improving knowledge of, and exciting applications of, the physics of materials at high density. 
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Los Alamos National Laboratory under Contract DE-AC52-06NA25396. This work is also 
supported by LDRDs 08-Em-002, 08-LW-004, and the National Laboratory User Facility 
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FIG. 1: (Color online) Compton scattering spectrum from laser compressed Fermi-degenerate 
beryllium. The width of the Compton feature is sensitive to density and its relative intensity is 
sensitive to the ion temperature. 
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FIG. 2: Phase diagram of hydrogen. Physical regimes are indicated with sohd hnes showing the 
plasma coupling parameter (F = 1) and the electron degeneracy parameter (9 = 1). For F > 178, 
the plasma freezes into a bcc Coulomb solid. The melting curve of H2 is also indicated for log T < 3. 
The dashed curve shows the dissociation and ionization boundaries in the low density gas. Above 
logP ~ 0.7, hydrogen is fully ionized. The Jupiter isentrope is shown by the heavy solid line. 
The regions probed by single, double and triple shock experiments on deuterium are indicated 
with dotted lines. Filled squares show the near isentropic compression data of Fortov et al. 
suggesting a PPT in hydrogen. The single and double shock helium points of Nellis et al. 3] are 
indicated with open squares. Finally, the hashed region outlines the locus of the shocked states 
achieved by Eggert et al. [13] from pre-compressed He targets. 
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FIG. 3: (Color online) Comparison of between theory (solid lines) and laser shock wave exper- 
iments (symbols). Helium was exposed to extreme temperatures and pressure that are relevant 
for planetary interiors. The colors represent different precompression ratios. The ability to pre- 
compress samples statically before launching the shock is an important experimental improvement 
that allows to probe deeper in the giant planet interiors. Good agreement between theory and 
experiment is found for the higher precompressions. 



16 




6300-6800 K 
2 Mbar 



\Rocks+/ 
13300 KVi^/l 5000-21 000 K 
38 Mbar 40 Mbar 

'0-10 Mr 



111 ]. Running along the left in 



FIG. 4: Schematic interior view of Jupiter, based on Militzer et al. 
black are pressures and temperatures from their model at three locations, as well as the core mass 
estimate {^15M^). The transition from molecular hydrogen (H2) to liquid metallic hydrogen (H^ 



is continuous. Running along the right in gray are these same estimates from Guillot j43l |. 
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FIG. 5: A comparison of theoretical mass-radius curves for gas giant planets and 40 observed 
transiting planets, using the models from Fortney et al. [2J]. The majority of these planets orbit 
at distances of only 0.02 to 0.05 AU from their parent stars, while the Earth orbits at 1 AU (by 
definition). The x-axis is mass in Jupiter masses, and the y-axis radius in Jupiter radii. The top 
two solid black curves are for pure H-He, 4.5 Gyr-old, giant planets at 0.02 AU and 0.045 AU 
from the Sun. (The Earth-Sun distance is 1 AU.) The thick dash-dot curve also shows models at 
0.045 AU, but with 25 (0.08 Mj) of heavy elements (ice+rock) in a core. A mass-radius curve 
for pure water planets is also shown. Gray diamonds are, left to right, Uranus, Neptune, Saturn, 
and Jupiter. Black diamonds with error bars are the transiting planets. Curves of constant bulk 
density (in g cm^'^) are overplotted in dotted gray. 
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FIG. 6: (Color online) Resolution tests obtained with 60 keV x-rays produced by a short pulse 
irradiated W-wire. 
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FIG. 7: (Color online) 60 keV radiograph of a shock compressed iron target. 
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FIG. 8: (Color online) Left: P-T profiles of terrestrial super-Earths. The family of planets with 
1-10 have a similar Fe/Si ratio as Earth. The highest internal P is 1.56 Mbar (156 GPa). The 
different phase transitions in the mantle are shown in dashed lines ranging from olivine (ol), wad- 
sleyite (wd) and ringwoodite (rw), perovskite (pv) and magnesiowusite (wu), and post-perovskite 
(ppv) and wu. The discontinuities are caused by the boundary layers at the top and bottom of 
the mantle. The mantles of super-Earths with masses larger than ~4 are mostly composed of 
ppv-|-wu, compared to the dominance of pv+wu on Earth. Right: Ternary Diagram for a 5 
planet. The radius of a planet with each mixture is shown in color with the color bar spanning the 
radius of the smallest (a 1 pure Fe planet with R = 5400 km) and largest (a 10 - pure 
H2O planet with R = 16000 km) super-Earth. The shaded region shows the unlikely compositions 
that can form a super-Earth from solar nebula condensation and secondary accretion constraints. 
A ternary diagram exists for every planetary mass value. 
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